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Abstract 

The  paper  presents  applications  of  the  large  eddy 
simulation  (LES)  methodology  on  the  Sandia/Sydney 
turbulent  blufjbody  burner  and  gas  turbine  combustors. 
LES  of  the  bluffbody  flame  is  performed  using  the  filtered 
density  function  (FDF)  submodel  and  a comprehensive 
augmented  chemical  mechanism  for  the  first  time.  The 
FDF  submodel  is  a sophisticated  turbulent-combustion 
submodel  that  directly  computes  the  joint  probability 
density  function  (PDF)  of  scalars  and  is  therefore 
considered  to  be  more  accurate  than  conventional 
assumed-PDF  type  models.  The  chemical  kinetics 
mechanism  involves  19  species  and  15  reaction-steps. 
The  mechanism  contains  both  Cl  and  C2  species  and  also 
involves  NO  formation  steps.  Owing  to  the  complexity  of 
the  mechanism,  numerical  integration  of  the  kinetics 
equations  is  performed  using  the  in  situ  adaptive 
tabulation  (ISAT)  scheme.  Mean  velocity  and 
species/temperature  fields  are  presented  and  compared  to 
experimental  data.  Results  show  that  the  computations 
are  in  good  agreement  with  data.  The  paper  also 
presents  LES  of  a gas  turbine  combustor.  LES  is 
performed  using  an  assumed  FDF  turbulent-combustion 
model  in  conjunction  with  the  fiamelet-generated 
manifold  method.  The  advantage  of  this  approach  is  that 
the  chemical  reaction  is  parameterized  by  only  two 
variables,  mixture  fraction  and  progress  variable.  Thus 
calculations  are  significantly  faster  than  those  with  the 
transport  FDF  model.  Circumferentially  averaged 
combustor  exit  fuel-air  ratio  profiles  are  compared  to 
measurement  data  for  two  liner  port  patterns.  It  is  shown 
that  the  LES  calculations  are  in  reasonable  agreement 
with  data  and  superior  to  Reynolds  averaged  Navier- 
Stokes  calculations.  These  calculations  indicate  that  LES 
of practical  combustion  systems  are  feasible  economically 
and  can  be  used for  design  analyses  more  routinely. 
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1.  Introduction 

Traditionally,  computational  modeling  of  turbulent 
reacting  flows  has  been  mainly  performed  using 
Reynolds-Averaged  Navier-Stokes  (RANS)  methods. 
However,  deficiencies  of  RANS  methods  for  swirling 
recirculating  flows  including  gas  turbine  combustor  flows 
are  becoming  increasingly  apparent.  LES  offers 
significant  advantages,  in  that,  it  captures  turbulent 
mixing  and  unsteady  flow  features  accurately.  In  LES,  a 
filtering  operation^’^  is  performed  on  the  governing 
equations  such  that  only  large-scale  flow  features  are 
resolved  on  the  computational  mesh,  while  the  small-scale 
features  are  modeled.  The  advantage  of  this  approach  is 
that,  if  one  could  afford  a fine-enough  mesh,  then 
theoretically,  very  little  modeling  would  be  required  in 
the  calculation.  The  disadvantage  of  this  method  is  that  a 
typical  LES  computational  mesh  is  significantly  larger 
than  a typical  RANS  mesh,  and  computationally 
expensive,  relatively.  However,  advances  in 
computational  modeling  and  the  availability  of  cheap 
distributed  computing  platforms  have  made  it  affordable 
to  use  LES  for  complex  geometries'^^“^'.  In  reeent  years, 
LES  has  become  a serious  candidate  for  analyzing 
reacting  flows  where  the  main  focus  of  interest  is  in 
regions  away  from  solid  boundaries. 

Although  LES  produces  better  predictions  of 
isothermal  flows  as  compared  to  RANS,  modeling  is  still 
needed  in  LES  for  accurate  predictions  of  combusting 
flows.  Combustion  occurs  at  the  very  small  scales  that 
are  not  resolved  on  the  computational  mesh.  Hence,  in 
addition  to  the  subgrid  turbulence  models,  subgrid 
combustion  models  are  also  required  to  account  for  the 
influence  of  unresolved  scales  of  the  flow  on  the  mean  (or 
filtered)  reaction  rates.  There  are  several  combustion 
models  proposed  for  LES.  Assumed  FDF  methods^^'  are 
the  most  widely  used  because  of  their  fast  turnaround 
time.  In  these  methods,  the  reaction  is  parameterized  by 
one  or  two  variables  such  as  mixture  fraction  and/or 
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progress  variably.  The  shape  of  the  PDF  of  these 
variables  is  assumed  a priori.  With  the  information 
available  on  the  PDF,  the  mean  reaction  rates  can  be 
easily  computed.  However,  applications  of  assumed-FDF 
methods  are  limited  to  reacting  flows  that  are  far  from 
extinction/ignition.  This  limitation  can,  in  principle,  be 
overcome  by  transport  PDF  methods^’“’l  These  methods 
are  considered  to  be  the  most  sophisticated  and  can  be 
applied  to  a variety  of  reacting  flows,  such  as  premixed, 
partially  premixed  and  nonpremixed  combustion  systems. 
In  transport  PDF  methods,  the  joint  PDF  of  reacting 
scalars  is  directly  computed  via  a PDF  transport  equation. 
The  advantage  of  these  methods  is  that  the  closure  term 
related  to  reaction  rates  is  obviated.  Also,  chemical 
kinetics  equations  of  arbitrary  complexity  can  be  easily 
used.  Due  to  the  complexity  of  performing  LES/FDF 
calculations,  applications  of  this  methodology  to  practical 
flames  have  been  limited.  Nevertheless,  there  is  a need  to 
fully  understand  the  predictive  capabilities  of  this  model. 

This  paper  presents  an  evaluation  of  the  PDF  model 
on  a standard  benchmark  bluffbody  stabilized  flame 
(Sandia/Sydney  Hydrogen/Methane  flame  [HMl]^'*'')  with 
a comprehensive  augmented  reduced  mechanism.  This 
mechanism  allows  the  computations  of  major  and  minor 
species  and  overall  provides  a more  realistic 
representation  of  the  chemical  processes  in  comparison  to 
simple  global  mechanisms.  However,  due  to  the  large 
number  of  species  and  reaction  steps  involved,  integration 
of  complex  mechanisms  is  computationally  expensive. 
The  in  situ  adaptive  tabulation  (IS AT)  scheme^'''  helps  to 
alleviate  this  problem  and  allows  the  use  of  such 
mechanisms  at  an  affordable  cost.  The  paper  presents  the 
velocities,  species  and  temperature  fields  fi'om 

calculations  and  compares  them  with  experimental  data. 

The  second  part  of  the  paper  presents  LES  of  a gas 
turbine  combustor.  The  subgrid  combustion  model  used 
in  these  calculations  is  based  on  an  assumed  PDF  method. 
LES  is  linked  to  a spray  transport  and  evaporation  model 
and  simulations  of  two  combustor  configurations  are 
presented.  Circumferentially  averaged  fuel-air  ratio 

profiles  at  the  combustor  exit  are  compared  with  rig 
measurements. 

The  paper  is  organized  as  follows.  Section  2 of  the 
paper  discusses  the  numerical  methodology  and 
combustion  models.  Section  3 presents  results  for  the 
bluffbody  flame.  Section  4 presents  LES  of  the  gas 
turbine  combustors  and  Section  5 presents  the  overall 
findings  of  this  study. 

2.  Numerical  Methodology  and  Combustion 
Models 

The  LES  flow  solver  is  based  on  a pressure-based, 
low  Mach  number,  collocated,  finite-volume  solver^^^,  to 


compute  velocity  fields.  The  numerical  scheme  is  an 
energy-conserving  scheme  that  is  second  order  accurate  in 
space  and  time.  Convective  and  diffusive  terms  in  the 
momentum  equations  are  discretized  using  a central 
scheme  and  time  derivative  is  discretized  using  a second- 
order  backward  difference  scheme.  The  subgrid  scale 
stress  model  used  is  based  on  the  dynamic  model. 
Convective  terms  in  the  scalar  transport  equations  are 
discretized  using  the  hybrid  linear-parabolic 
approximation  scheme  in  order  to  prevent  unphysical 
overshoots  and  undershoots  of  the  scalar  fields.  The  flow 
solver  is  linked  to  combustion  models  described  below. 


A.  FDF  Transport  Model 

The  FDF  transport  equation  is  given  as  References  7- 
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In  Eqs.  (1)  and  (2),  Fl  is  the  filtered  density  function, 
^1^  is  the  conditional  expectation,  ^ ^js  the  Favre 

filtered  mean,  p is  the  density,  S is  the  Dirac-delta 
function,  w,  is  the  instantaneous  velocity,  is  the 
diffusional  flux  of  species  a,  Sa  is  the  chemical  source 
term  for  species  a,  y/  is  the  chemical  composition  in 

scalar  space,  <j>  is  the  chemical  composition  in  physical 

space,  and  G is  the  filter  kernel.  In  Eq.  (1),  the  molecular 
mixing  term  (first  term  on  the  right  side)  is  typically 
modeled  using  the  interaction  by  exchange  with  mean 
mixing  model.  The  turbulent  convection  term  (second 
term  on  the  right  side)  is  modeled  by  the  gradient 
diffusion  model.  In  Eq.  (1),  it  can  be  seen  that  the  last 
term  on  the  right  side  is  due  to  chemical  reactions  and  is 
in  closed  form.  Thus  reactions  of  arbitrary  complexity 
can  be  implemented  to  describe  the  chemical  processes. 
Given  this  advantage  of  the  FDF  method  and  also  the  fact 
that  the  FDF  is  not  specified  a priori  as  in  assumed-FDF 
models,  it  is  expected  that  the  LES/FDF  method  will  yield 
better  predictions  of  not  only  overall  turbulent  flame 
structures,  but  also  local  extinction  and  reignition 
scenarios  in  these  flames. 

The  large  dimensionality  of  Eq.  (1)  precludes  the  use 
of  finite-difference  or  finite-volume  schemes  for  the 
solution  of  this  equation.  Monte-Carlo  methods  are 
typically  used  to  solve  such  equations.  In  Monte-Carlo 
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methods,  the  PDF  in  each  computational  cell  is 
represented  by  an  ensemble  of  notional  particles,  each 
with  a set  of  scalars  and  position  vector.  The  transport  of 
these  particles  in  physical  and  compositional  space 
mimics  the  solution  of  Eq.  (1).  The  mean  velocities 
needed  for  the  transport  of  these  particles  are  computed 
by  an  Eulerian  flow  solver.  Favre  filtered  scalar  values  in 
a cell  are  computed  by  ensemble  averaging  the  particle 
scalar  values  in  the  cell. 

B.  Assumed-FDF  Model  with  the  Flamelet- 
Generated  Manifold  Method 

In  contrast  to  transported  FDF  methods,  in  assumed 
FDF  methods^®',  the  shape  of  the  FDF  is  assumed  a priori. 
In  the  present  assumed  FDF  method,  two  scalar  variables 
are  considered — fuel  mixture  fraction  and  a progress 
variable.  This  allows  the  use  of  this  method  for  partially 
premixed  flames  commonly  observed  in  gas  turbine 
combustion  systems.  The  progress  variable  is  typically 
considered  to  be  a linear  combination  of  species.  In  the 
present  study,  the  progress  variable  is  considered  to  be 
CO2.  Within  the  flow-solver,  transport  equations  for 
mixture-fraction,  mixture-fraction  variance,  progress- 
variable  and  progress-variable  variance  are  solved.  The 
mixture  fraction  FDF  is  represented  by  a beta  function 
parameterized  by  mean  and  variance.  The  progress- 
variable  FDF  is  represented  by  a tri-delta  function,  also 
parameterized  by  mean  and  variance.  The  chemical 
source  terms  in  the  progress-variable  and  variance  of 
progress-variable  transport  equations  are  computed  via 
the  flamelet-generated  manifold  methodology^'^'.  In  this 
methodology,  it  is  assumed  that  a low-dimensional 
reaction  manifold  exists  in  the  composition  space,  and 
that  a reaction  path  would  move  towards  the  manifold  and 
subsequently  remain  on  this  manifold.  The  manifold  is 
generated  by  performing  unstrained  premixed  flame 
calculations  for  different  mixture-fraction  values  using 
detailed  chemical  kinetics.  The  chemical  structure  of  the 
flame  is  parameterized  by  mixture  fraction  and  progress 
variable.  The  advantage  of  this  approach  is  that  the 
source  terms  can  be  pre-tabulated  (before  the  LES  is 
performed)  in  terms  of  mean  and  variance  of  mixture- 
fraction  and  progress-variable.  During  the  LES 
calculation,  linear  interpolation  can  be  used  to  compute 
the  source  terms.  The  assumed  FDF  method  is  therefore 
computationally  cheap  and  well  suited  for  routine 
practical  applications. 


3.  Sandia/Sydney  Bluff-Body  Stabilized 
Flame  (HMl) 

The  HMl  burner''®'  consists  of  a central  jet  of 
methane/hydrogen  (1;1  volume  ratio)  surrounded  by  a 
bluffbody.  The  bluffbody  is  immersed  in  a coflowing  air 
stream.  The  diameter  of  the  central  jet  and  the  bluffbody 
are  3.6  mm  and  50  mm,  respectively.  The  bulk  velocities 
of  the  two  streams  are  1 18  m/s  and  40  m/s,  respectively. 

LES  of  the  burner  is  performed  on  a 136  x 102  x 102 
computational  mesh  (1.4  million  grid  points)  that  starts  at 
the  face  of  the  bluffbody.  A fully  developed  pipe  flow 
profile  is  imposed  at  the  fuel  stream  boundary  in  order  to 
closely  capture  the  penetration  of  the  fuel  jet.  At  all  the 
air  stream  boundaries,  a plug  flow  profile  is  imposed. 
Further  studies  are  needed  to  understand  the  sensitivity  of 
the  predictions  on  the  profiles  imposed  on  the  fuel  and  air 
streams.  The  FDF  transport  model  in  conjunction  with  an 
augmented  reduced  mechanism  for  methane  is  used  in  the 
calculation  of  this  flame.  In  order  to  balance 
computational  effort  and  accuracy,  five  particles  per  cell 
are  used  in  the  simulation.  However,  to  minimize 
stochastic  errors,  all  particles  in  neighboring  cells  are  also 
considered  in  the  computation  of  ensemble  means. 

The  augmented  scheme  (Table  1)  is  based  on  a 19- 
species  15-step  Cl  and  C2  mechanism''^'.  Due  to  the 
complexity  of  this  mechanism,  integration  of  the  kinetics 
equations  is  performed  using  ISAT.  The  objective  of  the 
ISAT  scheme  is  to  efficiently  integrate  the  coupled  set  of 
reaction  rate  equations: 

^^  = 5Jd)(0]  (3) 

i.e.,  given  the  initial  condition  <t(/o)>  the  reaction  rate 
equation  is  integrated  for  a time  O(to  + A/).  A table  is 
built  as  the  reactive  flow  calculation  progresses,  so  that 
only  the  accessed  region  of  the  composition  space  is 
tabulated.  Since  the  accessed  region  is  a small  subset  of 
the  realizable  region,  the  dynamic  creation  of  the  table 
enables  more  efficient  integration  of  the  reaction 
equations  than  when  the  table  is  created  a priori  and  the 
entire  realizable  region  is  tabulated.  The  accuracy  of  the 
ISAT  integration  is  determined  by  a user  specified 
tolerance  parameter  e.  In  the  present  computations,  a 
value  of  10"“'  was  chosen  for  s,  to  achieve  acceptable 
accuracy  levels.  Calculations  were  performed  on  nine 
processors  on  a Linux  machine  and  require  approximately 
a month  for  full  convergence  of  statistical  quantities.  The 
majority  of  the  time  was  spent  on  chemistry  computations 
and  in  achieving  convergence  of  the  statistical  averages 
including  those  of  trace  species  such  as  CO  and  NO. 
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Table  1. 19-Species,  15-Step  Augmented  Reduced 
Mechanism 

H+0.502=OH 

H2+0.502=H+OH 

H02=0.502+OH 

0.502+H202=OH+H02 

0.502+0.5C2H2=H+CO 

CH3+CO+C2H4=0.502+CH4+1  .5C2H2 

0.502+2CH3=H2+CH4+CO 

0.502+CH3=H+CH20 

0.502+CH4=OH+CH3 

0.502+CO=C02 

0.502+C2H6=CH4+CH20 

H+0H=H20 

H+CH4+NO+HCN=0.502+2CH3+N2 
H+0.502+CH4+HCN=2CH3+NO 
0.502+CH4+NH3+HCN=H20+2CH3+N2 


Figure  1 presents  the  instantaneous  and  time- 
averaged  snapshot  of  the  computed  temperature  field. 
The  instantaneous  temperature  field  shows  an  intermittent 
flame  structure.  The  time-averaged  field  shows  high 
temperatures  in  the  near-field  region  and  low  temperature 
downstream.  The  near  field  of  the  bluflbody  flame  is 
characterized  by  two  counter-rotating  vortex  rings  (not 
shown  in  the  figure)  that  help  to  stabilize  the  flame. 
Figure  2 shows  the  axial  and  radial  velocities, 
respectively,  at  several  axial  stations  in  the  flowfield. 
Note  that  the  velocities  are  reasonably  well  predicted. 
Some  discrepancies  are  seen  in  the  radial  velocities  at  far 
downstream  locations  that  may  be  related  to  the 
coarseness  of  the  grid  near  the  exit  boundary.  Figure  3 
shows  the  radial  profiles  of  temperatures.  In  the  near- 
field region,  the  temperatures  are  well  represented; 
however,  at  x = 120  mm,  peak  temperatures  are 
underpredicted.  Similar  features  are  also  seen  in  the 
radial  profiles  of  H2,  CO,  CO2  and  NO  mass  fi-actions 
(Figure  4).  Overall,  it  can  be  seen  that  the  predictions 
reproduce  the  data  well.  The  near-field  region  is  well- 
represented;  however,  some  discrepancies  are  seen  in  the 
far-field  region,  which  may  be  mitigated  by  increasing  the 
grid  resolution. 

4.  Gas  Turbine  Combustor 

The  application  of  LES  to  practical  combustion 
systems  is  demonstrated  by  performing  calculations  of  a 
gas  turbine  combustor.  The  computational  domain 
consists  of  the  region  within  the  combustor  liner  only. 
The  domain  starts  downstream  of  the  swirler  vanes.  The 
swirl  boundary  conditions  and  liner  primary  and  dilution 
port  boundary  conditions  were  determined  by  network 
analyses  using  measured  pressure  measurements  on  the 
combustor  test-rig.  The  grid  used  in  the  calculations  has 


approximately  3 million  hexahedral  elements  distributed 
across  20  processors  on  a Linux  machine.  The  overall 
duration  for  each  calculation  was  about  three  weeks. 

The  flow  solver  was  linked  to  the  assumed  FDF 
combustion  model  described  in  Section  2 (B).  Since 
liquid  fuel  is  used  in  the  combustor,  a Lagrangian  spray 
droplet  transport  and  evaporation  modeF'  was  linked  to 
the  flow  solver.  Discrete  size  droplets  were  injected 
downstream  of  the  fiiel  nozzle.  These  droplets  were 
transported  by  the  surrounding  liquid  and  evaporated  by 
the  ambient  air  stream.  The  transport  equations  reflected 
the  impact  of  droplet  evaporation  via  source  terms. 

Two  different  liner  wall  port  patterns  (Mod  A and 
Mod  B)  were  analyzed  (Figure  5).  Mod  A had  two  rows 
of  holes  on  both  inner  and  outer  liner  walls.  Mod  B had 
three  rows  of  holes  on  the  inner  wall.  Also,  the  size  of  the 
holes  in  the  second  row  of  the  outer  liner  wall  was 
reduced  fi’om  Mod  A.  Figures  6 and  7 present  the 
instantaneous  contours  of  velocity  magnitude  and 
temperature  in  Mod  A and  Mod  B,  respectively.  The 
dynamical  features  of  the  flow  and  reacting  fields  can  be 
clearly  seen.  In  comparison  to  Mod  A,  the  higher  airflow 
from  the  inner  wall  in  Mod  B results  in  lower 
temperatures  close  to  the  inner  wall.  Similarly,  the  lower 
airflow  from  the  outer  wall  in  Mod  B results  in  higher 
temperatures  close  to  the  outer  wall.  Figure  8 presents 
circumferentially  averaged  fuel/air  ratio  versus  span  at  the 
combustor  exit  from  measurements  and  computations.  It 
can  be  seen  that  the  LES  calculations  reproduce  the 
experimental  trend  reasonably  well.  In  Mod  A,  RANS 
was  also  used  to  compute  the  exit  profile.  It  can  be  seen 
from  the  figure  that  the  prediction  from  LES  is  superior  to 
that  from  RANS. 

5.  Summary 

The  paper  presents  applications  of  LES  to  laboratory 
flames  and  practical  gas  turbine  combustors.  LES  of  the 
bluffbody-stabilized  turbulent  hydrogen-methane  flame 
was  performed  using  the  transport  FDF  model  with  a 
comprehensive  augmented  reduced  mechanism  that 
involved  19  species  and  15  reaction-steps  for  the  first 
time.  The  chemical  mechanism  involved  both  Cl  and  C2 
species  and  was  integrated  using  ISAT.  Mean  velocity 
and  species/temperature  fields  were  presented  and 
compared  to  experimental  data.  Results  showed  that  the 
computations  are  in  good  agreement  with  data.  The  paper 
also  presents  LES  of  a gas  turbine  combustor. 
Simulations  of  the  combustor  were  performed  using  an 
assumed  FDF  turbulent-combustion  model  in  conjunction 
with  the  flamelet-generated  manifold  method. 
Circumferentially  averaged  combustor  exit  fuel-air  ratio 
profiles  were  compared  to  measurement  data  for  two  liner 
port  patterns.  It  was  shown  that  the  LES  calculations 
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were  in  reasonable  agreement  with  data.  These 
calculations  indicated  that  LES  of  practical  combustion 
systems  is  feasible  and  can  be  routinely  used  for  design 
analyses. 
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Figure  1.  Computed  instantaneous  (left)  and  time- 
averaged  (right)  temperature  fields  in  the  HM1  burner 
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Figure  2.  Comparison  of  axial  (top)  and  radial  (bottom) 
velocities  with  data  in  the  HM1  burner 
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Figure  3.  Comparison  of  predicted  temperature 
against  data  for  the  HM1  burner 


Figure  4.  Comparison  of  predicted  H2,  CO,  CO2  and  NO 
mass  fractions  against  data  for  the  HM1  burner 
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Figure  5.  The  two  combustor  liner  hole  configurations 
investigated;  Mod  A (left),  Mod  B (right) 


Figure  6.  LES  of  Mod  A combustor:  instantaneous 
velocity  magnitude  (top  left);  instantaneous 
temperature  (top  right);  instantaneous  temperature  at 
combustor  exit  (bottom).  All  units  are  normalized. 


Figure  7.  LES  of  Mod  B combustor:  instantaneous 
velocity  magnitude  (top  left);  instantaneous 
temperature  (top  right);  instantaneous  temperature  at 
combustor  exit  (bottom).  All  units  are  normalized. 
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Figiree  8.  Circumferentially  averaged  normalized  fuel- 
air  ratio  versus  normalized  span  at  the  exit  of  the 
combustor:  Mod  A (left),  Mod  B (right) 
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